Internal Spatiotemporal Stochastic Resonance in a Microscopic Surface Reaction 

Model 

O. Kortliike* *, V.N. Kuzovkov'''*, and W. von Niessen* 
^ Institut fiir Physikalische und Theoretische Chemie, Technische Universitdt Braunschweig, Hans-Sommer-Strafle 10, 38106 

Braunschweig, Germany 
''Institute of Solid State Physics, University of Latvia, 8 Kengaraga Street, LV - 1063 RIGA, Latvia 

(Received February 1, 2008) 



in 



> 

in 

O 
0\ 
0^ 
-(— • 



o 



% 



PACS numbers: 05.40.Ca, 05.10.Ln, 05.45.-a, 05.45.Xt 

We show the existence of internal stochastic resonance in a 
microscopic stochastic model for the oscillating CO oxidation 
on single crystal surfaces. This stochastic resonance arises di- 
rectly from the elementary reaction steps of the system with- 
out any external input. The lattice gas model is investigated 
by means of Monte Carlo simulations. It shows oscillation 
phenomena and mesoscopic pattern formation. Stochastic 
resonance arises once homogeneous nucleation in the individ- 
ual surface phases (reconstructed and non-reconstructed) is 
added. This nucleation is modelled as a noise process. As a 
result, synchronization of the kinetic oscillations is obtained. 
Internal stochastic resonance may thus be an internal regula- 
tion mechanism of extreme adaptability. 

The term stochastic resonance (SR) is given to the 
somewhat counter-intuitive phenomenon that in a non- 
Unear system a weak signal can be amplified by the assis- 
tance of noise. It has been introduced in 1981 by Benzi 
et al. ||l|,|| in the context of a study about the period- 
ically recurrent ice ages. Over the last two decades it 
has continuously attracted increasing attention and was 
shown to occur in many syste ms p-l^ in biology |£|-||, 
chemistry |9|-[ll|, and physics |l2|-|l4[|. Generally sys- 
tems showing SR are described in a formal mathemat- 
ical way using phenomenological macroscopic equations 
of the mean field type including (i) a bistable system 
with an activation barrier or some sort of threshold, (ii) 
a weak coherent input, and (iii) a strong external noise 
which helps to overcome the activation barrier. These 
macroscopic equations in a sense are able to describe 
many different systems (because stochastic resonance is 
a general phenomenon which occurs in many natural sys- 
tems) and have been used in the description of stochas- 
tic resonance phenomena in Nd-YAG lase rs p5| | homo- 
geneous [16| as well as heterogeneou s p7 ,lS| chemical 
reactions, bistable quantum systems |19|] or the Lotka- 
Volterra model p^. More complex systems (e.g. a sum- 
ming network of excitable units pl| , sheep populations 
|P2| , two-dimensional excitable media showing spatio- 
temporal pattern formation [p3|, sensory systems in cray- 
fish [gj or in the visual cortex |p4| , or neuron- like systems 
Ipsi ) are generally modeled via Langevin equations or the 
Fitzhugh-Nagumo model. 

In addition a few special systems have been investi- 
gated in a more general manner via a macroscopic math- 
ematical description, e.g. an autonomous oscillating sys- 



tem p6J, a system in the limit of weak noise p^, a 
system showing stochastic multiresonance p8[ , and non- 
dynamical systems with both internal and external noise 
p9[ . Note that in the latter case the internal noise is 
modeled in the same way as an external noise and that it 
is only regarded as a general internal noise without spec- 
ifying the physical background. Computer simulations 
performed to date consider coupled neurons or general 
threshold devices, which are mesoscopic models, i.e. the 
microscopic physical picture is again neglected. 

Our present model system is very unusual in the re- 
search on SR. It gives for the first time a microscopic 
description of the phenomenon of internal stochastic res- 
onance and demonstrates its physical reasons on the mi- 
croscopic (atomic) length scale. It is based on stochastic 
transitions, each with a clear physical meaning. With- 
out noise, the system exhibits a spatially extended het- 
erogeneous stable state with inherent local oscillations. 
More important, the noise is not an external input but 
corresponds to a physically realistic internal nucleation 
process. Because of the clear physical picture on the 
microscopic (atomic) level our model is of course spe- 
cialized and cannot describe a large variety of different 
systems. But on the other hand the results of this model 
and the conclusions which can be drawn are very general 
ones and suggest that internal SR may be the reason for 
inherent synchronization and cooperative phenomena in 
many physical, chemical and biological systems. 

We consider a slightly modified version of a previ- 
ously presented model for the catalytic CO-l-1/2 O2 re- 
action on Pt single crystal surfaces |3y,^^, which shows 
different types of kinetic oscillations in agreement with 
experimental results. The model involves CO adsorp- 
tion, desorption and diffusion, dissociative O2 adsorp- 
tion and two surface phases (reconstructed and non- 
reconstructed) which form and propagate governed by 
the coverage with CO. The details are given below. An 
extended version of the model for the CO-t-NO reaction 
on Pt(lOO) is able to describe the experimentally ob- 
served transition into chaotical behavior via the Feigen- 
baum route [|2|. Our model follows the well known 
model by Ziff, Gulari and Barshad ||3^ (ZGB model) and 
is investigated by means of Monte Carlo (MC) simula- 
tions. The Pt(llO) surface of the catalyst is represented 
by a square lattice of side length L and lattice constant 
a = 1. From experiment p4| it is well known that kinetic 



oscillations are closely connected with the a ^ P recon- 
struction of the Pt(llO) surface, where a and /3 denote 
the 1x2 and the 1x1 surface phase, respectively. In 
our model CO is able to adsorb onto a free surface site 
with rate yco = V and to desorb from the surface with 
rate fc, independent of the surface phase the site belongs 
to. O2 adsorbs dissociatively onto two nearest neighbor 
(NN) sites with different sticking coefficients s onto the 
two phases {sa — 0.5, sp = 1). Therefore we get the 
oxygen adsorption rates TJo — ^ — y and j/q = 2(1 — y) 
for the a and (5 phase, respectively. For O2 adsorption 
directly at the phase border where one site belongs to the 
a and the other one to the (3 phase the geometric mean 
of these adsorption rates is used. In addition, CO is able 
to diffuse with rate D via hopping onto a vacant NN site. 
The CO+O reaction occurs, if CO hops to a site which 
is covered by O and the reaction product CO2 desorbs 
immediately from the surface. All these processes are 
associated with the above kinetic transition rates of the 
stochastic model which therefore determine the relative 
speed of the individual reaction processes. In the present 
study we use y = 0.51, D = 100, and k = 0.1 as standard 
values because CO diffusion is by far the fastest process. 
For details see refs. [|0|-§2|j35) . 

The a ^ /3 phase transition is modeled as a linear 
phase border propagation. Consider two NN surface sites 
in the state a/3. The transition a(3 — > aa {aP — > /3/3) oc- 
curs if none (at least one) of these two sites is occupied by 
CO. This phase border propagation mechanism mimicks 
the growth of the /3 phase because of the larger adsorption 
energy of CO on the 1x1 phase than on the 1x2 phase 
[Q . The individual phases are stable or metastable. The 
direct transition from a globally homogeneous a phase 
into a homogeneous /? phase (or vice versa) is impossi- 
ble; the activation barrier is infinite. The stability of the 
individual phases depends on the chemical coverage Qi of 
species i on the surface of the catalyst. For 8co < 0.3 the 
a phase, for larger values the p phase is stable |gl[] . The 
coverages of CO and O vary in the course of the reaction 
because of the different sticking coefficients of O2 on the 
two surface phases. Starting with a heterogeneous distri- 
bution of the a and /3 phase the activation barrier for the 
surface phase transition is finite, but the transition into 
a globally homogeneous phase does not occur because of 
the finite surface phase propagation velocity, i.e. the a or 
/3 phases cannot grow to macroscopic islands. Therefore 
the oscillations remain local, interfere and cancel each 
other on sufficiently large surfaces |^,|l[. The system 
exists in a heterogeneous, dynamically stable state with 
oscillations, which are locally synchronized by CO diffu- 
sion but disappear on the macroscopic length scale for 
large lattices. The a and /3 phase, however, build almost 
homogeneous islands on a mesoscopic length scale. 

The nucleation is modeled as a spontaneous a ^ (3 or 
/3 — > a transition of a single site, completely independent 



of its neighbors or the particle adsorbed onto this site. 
Therefore this nucleation is a homogeneous process which 
corresponds to a weak noise which generates dynamic de- 
fects on the surface. The term weak noise is used because 
the nucleation rate 7 S [10~^, 10~^] is for the relevant val- 
ues several orders of magnitude smaller than the other 
transition rates which are of the order of 10~^ to 10^. 
The defects grow or vanish via the a ^ f3 reconstruction 
depending on their chemical neighborhood, i.e. the pres- 
ence or absence of CO. It has to be emphasized that all 
phenomena such as local oscillations, growth and decline 
of the heterogeneously distributed surface phase islands, 
quasiperiodical and chaotical behavior exist even with- 
out the consideration of the nucleation process [p0|-p2[ . 
Because of the very small nucleation rate, which intro- 
duces only very few defects into the existing heteroge- 
neous surface state, one might suppose that the defects 
could have no influence on the oscillating system, but we 
shall demonstrate that the nucleation process may have 
a profound influence on the system behavior leading to 
cooperative phenomena and the synchronization of the 
local oscillations. 

In comparison with the standard problem of SR there 
exist several major differences, (i) We have a two-phase 
system with an additional chemical coverage instead of 
a bi-stable system. In addition none of the two phases 
becomes homogeneous. Both coexist in a dynamically 
stable heterogeneous state, (ii) No external coherent in- 
put is considered in our model. The local oscillations 
originate from the kinetic definition of the model itself 
and are very small on a macroscopic, global length scale. 
These internal oscillations can be seen as a substitute 
of the commonly used coherent input. The macroscopic 
synchronization of the local oscillations then corresponds 
to a large output signal, (iii) The source of noise is an 
internal physical process of the system, the homogeneous 
nucleation, which breaks the mesoscopic homogeneity of 
the surface phase islands. The noise process is therefore 
not an addition to a periodic input but is independent 
of the oscillations. In most previous studies a weak co- 
herent input is coupled with a strong external noise, (iv) 
The model is investigated via MC simulations which cor- 
respond most closely to a hierarchy of master equations 
with all correlations included and then mapped onto a 
finite lattice. The simulation procedure contains addi- 
tional noise by its very nature but such a noise does not 
lead to SR. We may characterize the structure of our 
model by saying that practically all processes which give 
rise to SR (except for the particle flux to the surface) are 
internal to the system. 

As can be seen in fig. P the nucleation (noise) has a 
strong influence on the system behavior. Without or 
with small nucleation rates 7 < 5 • 10~^ only local os- 
cillations exist. The amplitude of the global oscillations 
vanishes for simulations on large lattices pHmt]. With 
increasing nucleation rate (strength of the noise) the lo- 



cal oscillations are synchronized on a macroscopic scale 
and almost reach the theoretical maximum of S{uj) = 0.5 
(0/3 varies between and 1). This holds also for larger 
lattices up to L = 4096 which is the current limit for our 
simulations. Further increase of the noise decreases this 
synchronization until at a nucleation rate of 7 = 10^^ 
the system is completely governed by strong noise. This 
behavior is the fingerprint of SR 0] . 

In addition the noise has an influence on the frequency 
of the oscillations (see fig. 0) . The normal frequency of 
the system is coq, which can only be observed in simula- 
tions on small lattices. For nucleation rates 7 > 5 • 10^'' 
the phenomenon of SR occurs. The nucleation forces the 
system to oscillate with a different frequency uj starting 
at a; w 2/3 wq for j — 5 ■ 10~^. With increasing nucle- 
ation rate the frequency increases as well up to a value 
of about oj « 2^0 at 7 = 5 • 10~^. This increase in the 
frequency is based on the increasing number of dynamic 
phase defects which grow very fast and accelerate the cor- 
responding phase transition. If the number or density of 
defects becomes too large the oscillating behavior breaks 
down. The shape of the curve in fig. ^ in the interval 
7 S [10~^, 10~^] is in very good agreement with the one 
obtained in the study of a general autonomous system by 
Haken et al. ||||. 

The underlying mechanism of this internal stochas- 
tic resonance effect is unexpectedly simple. During the 
growth and decline of the individual surface phases only 
a few very small residual phase islands remain in do- 
mains where the local synchronization leads to large am- 
plitudes in the phase oscillations. These residual islands 
are spatially separated at a mean distance R^. Adsor- 
bate diffusion is well known to synchronize individual sur- 
face domains within the socalled synchronization lentgh 
^ oc \fljf fsy,^ , where D is the diffusion rate and T is 
the time for one oscillation period. Without nucleation 
^ < i?r holds and only locally synchronized oscillations 
exist. In addition, the residual islands can only com- 
bine into a homogeneous phase if the second condition 
i?r ~ VT for the phase border propagation is fulfilled. 
But this condition is also violated because R^ > VT 
holds. Nucleation and subsequent growth now leads to 
new small phase islands (see fig. H). The mean distance 
between the individual phase islands decreases to R^i, for 
which Rn < £, < Ri and Rn ~ VT holds for proper nu- 
cleation rates. Therefore the separated phase islands can 
now be connected via island growth and synchronized via 
CO diffusion. This results in macroscopic synchronized 
oscillations. 

The existence of the surface phase nucleation is well 
known but has not been investigated experimentally yet, 
in contrast to the initial growth of small surface phase 
domains p6| . The nucleation rate should depend on the 
temperature 7 — 7(r), but it is almost impossible to 
achieve an isolated variation of the nucleation rate un- 
der experimental conditions because all other parameter 



such as CO desorption and CO diffusion also strongly 
depend on the temperature. It might thus be very diffi- 
cult to experimentally verify the mechanism behind the 
SR phenomenon in our model but there should be other 
systems where it is feasible. 

The nucleation of dynamic surface defects as an inter- 
nal process generates globally synchronized oscillations 
in our CO-|-O2/Pt(110) model reaction system via SR. 
Moreover, because noise is always present in real sys- 
tems, this type of internal SR should be a very general 
phenomenon and it may be the reason for cooperative 
phenomena and internal synchronization via noise pro- 
cesses in many physical, chemical, and biological systems 
where is has not been investigated experimentally or even 
suspected to be present to date. This especially holds for 
systems which exhibit inherent oscillations which are syn- 
chronized on macroscopic length scales. In this case often 
noise effects are supposed to be negligible, but as shown 
above they can also be the origin of those cooperative 
phenomena. We believe that SR may thus be an inter- 
nal regulation mechanism of extreme adaptability. This 
conclusion is drawn to search for internal SR in a variety 
of experiments, because to date experiments designed to 
study the role of internal noise have been inconclusive 
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FIG. 1. Global amplitude S{uj) of the /3 phase coverage 
as a function of the nucleation rate 7 which gives the noise 
strength. A stochastic resonance maximum can be seen. The 
scatter of the points around the maximum are a result of the 
finite frequency interval due to FFT analysis. The values of 
the amplitude are averaged over 20 simulation results. All 
other parameters are kept constant nt y — 0.51 (CO adsorp- 
tion), k = 0.1 (CO desorption), D = 100 (CO diffusion), 
V = 1 (surface phase propagation), and L = 256. 




FIG. 2. Ratio of the frequencies uj/u>o as a function of 
the nucleation rate 7 which gives the noise strength. Each 
value is averaged over 20 independent MC simulations. The 
errors are smaller than the symbol size. All other parameters 
are kept constant at y — 0.51 (CO adsorption), k — 0.1 (CO 
desorption), D — 100 (CO diffusion), V = 1 (surface phase 
propagation), and L — 256. 
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FIG. 3. New phase nuclei lead to a decreased mean dis- 
tance _Rn between the individual phase islands. This then 
leads to global synchronization over the whole lattice. See 
text for details. 



